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Abstract. A brief review of the nanoscale free-electron model of metal nanowires is 
presented. This continuum description of metal nanostructures allows for a unified 
treatment of cohesive and conducting properties. Conductance channels act as 
delocalized chemical bonds whose breaking is responsible for jumps in the conductance 
and force oscillations. It is argued that surface and quantum-size effects are the two 
dominant factors in the energetics of a nanowire, and much of the phenomenology of 
nanowire stability and structural dynamics can be understood based on the interplay of 
these two competing factors. A linear stability analysis reveals a sequence of "magic" 
conductance values for which the underlying nanowire geometry is exceptionally 
stable. The stable configurations include Jahn- Teller deformed wires of broken axial 
symmetry. The model naturally explains the experimentally observed shell and 
supershell structures. 
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1. Introduction 

The past decades have seen an accelerating miniaturization of both mechanical 
and electrical devices, so that a better understanding of properties of ultrasmall 
systems is required in increasing detail. The first measurements of conductance 
quantization in the late 1980s ( IWees et al. 19881 IWharam et al. 19881) in constrictions 



of two-dimensional electron gases formed by means of gates, have demonstrated the 
importance of quantum confinement effects in these systems and opened a wide field 
of research. A major step has been the discovery of conductance quantization in 
metallic nanocontacts ( Agrait et al. 1993 Brandbyge et al. 1995 IKrans et al. 19951) : 



The conductance measured during the elongation of a metal nanowire is a steplike 
function where the typical step height is frequently near a multiple of the 
conductance quantum Go = 2e^ jh^ where e is the electron charge and h Planck's 
constant. Surprisingly, this was initially not interpreted as a quantum effect but 
rather as a consequence of abrupt atomic rearrangements and elastic deformation 
stages. This interpretation, supported by a series of molecular dynamics simulations 
f lLandman et al. 19901 ITodorov et al. 19981) . was claimed to be confirmed by another 
pioneering experiment ( IRubio et al. 19961 IStalder et al. 19961) measuring simultaneously 



the conductance and the cohesive force of gold nanowires with diameters ranging from 
several Angstroms to several nanometers. As the contact was pulled apart, oscillations 
in the force of order InN were observed in perfect correlation with the conductance 
steps. 

It came as a surprise when Stafford et al. (1997) introduced the free-electron model 
of a nanocontact - referred to as the Nanoscale Free- Electron Model (NFEM) henceforth 
- and showed that this comparatively simple model, which emphasizes the quantum 
confinement effects of the metallic electrons, is able to reproduce quantitatively the main 
features of the experimental observations. In this approach, the nanowire is understood 
to act as a quantum waveguide for the conduction electrons (which are responsible for 
both conduction and cohesion in simple metals): Each quantized mode transmitted 
through the contact contributes Gq to the conductance and a force of order Ep/Xp to 
the cohesion, where Ep and Xp are the Fermi energy and wavelength, respectively. 
Conductance channels act as delocalized bonds whose stretching and breaking is 
responsible for the observed force oscillations, thus explaining straightforwardly their 
correlations with the conductance steps. 

Since then, free-standing metal nanowires, suspended from electrical contacts at 
their ends, have been fabricated by a number of different techniques. Metal wires 
down to a single atom thick were extruded using a scanning-tunneling microscope 
tip fIRubio et al. 19961 lUntiedt et al. 1997p . Metal nanobrid ges were shown to "self- 



assemble" under electron-beam irradiation of thin metal films (Kondo et al. 1997, 
2000, Rodrigues et al. 2000), leading to nearly perfect cylinders down to four atoms 
in diameter, with lengths up to fifteen nanometers. In particular, the mechanically- 
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controllable break junction technique, introduced by Moreland and Ekin (1985) and 
refined by Ruitenbeek and coworkers (IMuUer et al. 19921) . has allowed for systematic 
studies of nanowire properties for a variety of materials. For a survey see the review by 
Agrait et al. (2003). 

A remarkable feature of metal nanowires is that they are stable at all. Most 
atoms in such a thin wire are at the surface, with small coordination numbers, so 
that surface effects play a key role in their energetics. Indeed, macroscopic arguments 
comparing the surface- induced stress to the yield strength indicate a minimum radius 
for solidity of order ten nanometers (Zhang et al. 2003). Below this critical radius 
and absent some other stabilizing mechanism, plastic fiow would lead to a Rayleigh 
instability ( IChandrasekhar 198ip breaking the wire apart into clusters. Already in 
the 19th century Plateau (1873) realized that this surface-tension-driven instability is 
unavoidable if cohesion is due solely to classical pairwise interactions between atoms. 
The experimental evidence accumulated over the past decade on the remarkable stability 
of nanowires considerably thinner than the above estimate clearly shows that electronic 
effects emphasized by the NFEM dominate over atomistic effects for sufficiently small 
radii. 

A series of experiments on alkali metal nanocontacts (Yanson et al. 1999, 2001) 
identified electron- shell effects, which represent the semiclassical limit of the quantum- 
size effects discussed above, as a key mechanism influencing nanowire stability. 
Energetically-favorable structures were revealed as peaks in conductance histograms, 
periodic in the nanowire radius, analogous to the electron-shell structure previously 
observed in metal clusters (Ide Heer 1993p . A supershell structure was also observed 
( lYanson et al. 20"00ll . in the form of a periodic modulation of the peak heights. More 
recently, such electron-shell effects have also been observed, even at room temperature, 
for the noble metals gold, copper, and silver (Diaz et al. 2003, Mares et al. 2004, 2005) 
as well as for aluminum ( IMares et al. 200"7ll . 

Soon after the flrst experimental evidence for electron shell effects in metal 
nanowires, a theoretical analysis using the NFEM found that nanowire stability can 
be explained by a competition of the two key factors, surface tension and electron-shell 
effects (IKassubek et al. 200ip . Both linear ( |Zhang et al. 20031 lUrban et al. 2003!) and 
nonlinear (Biirki et al. 2003, 2005a) stability analyses of axially symmetric nanowires 
found that the surface-tension driven instability can be completely suppressed in the 
vicinity of certain "magic radii." However, the restriction to axial symmetry implies 
characteristic gaps in the sequence of stable nanowires, which is not fully consistent 
with the experimentally observed nearly perfect periodicity of the conductance peak 
positions. A Jahn- Teller deformation breaking the symmetry can lead to more stable 
deformed conflgurations. Recently, the linear stability analysis was extended to wires 
with arbitrary cross-section (Urban et al. 2004a, 2006). This general analysis conflrms 
the existence of a sequence of magic cylindrical wires of exceptional stability which 
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represent roughly 75% of the main structures observed in conductance histograms. The 
remaining 25% are deformed and predominantly of elliptical or quadrupolar shapes. This 
result allows for a consistent interpretation of experimental conductance histograms for 
alkali and noble metals, including both the electronic shell and supershell structures 
flUrban et al. 2nn4bl) . 

This chapter is intended to give an introduction to the NFEM. Section[2]summarizes 
the assumptions and features of the model while the general formalism is described in 
Sec. [3l In the following sections, two applications of the NFEM will be discussed: First, 
we give a unified explanation of electrical transport and cohesion in metal nanocontacts 
(Sec. Hj) and second, the linear stability analysis for straight metal nanowires will be 
presented (Sec. [5]). The latter will include cylindrical wires as well as wires with broken 
axial symmetry, thereby discussing the Jahn- Teller-effect. 



2. Assumptions and limitations of the NFEM 

Guided by the importance of conduction electrons in the cohesion of metals, and by the 
success of the jellium model in describing metal clusters ( Ide Heer 1993t IBrack 1993p . 



the NFEM replaces the metal ions by a uniform, positively charged background that 
provides a confining potential for the electrons. The electron motion is free along the 
wire, and confined in the transverse directions. Usually an infinite confinement potential 
(hard-wall boundary conditions) for the electrons is chosen. This is motivated by the 
fact that the effective potential confining the electrons to the wire will be short ranged 
due to the strong screening in good metals. 

In a first approximation electron-electron interactions are neglected, which is 
reasonable due to the excellent screening (IKassubek et al. 19991) in metal wires with 



G > Gq. It is known from cluster physics that a free electron model gives qualitative 
agreement and certainly describes the essential physics involved. Interaction, exchange 
and correlation effects as well as a realistic confinement potential have to be taken into 
account, however, for quantitative agreement]^ From this we infer that the same is 
true for metal nanowires, where similar confinement effects are important. Remarkably, 
the electron-shell effects crucial to the stabilization of long wires are described with 
quantitative accuracy by the simple free-electron model, as discussed below. 

In addition, the NFEM assumes that the positive background behaves like an 
incompressible fluid when deforming the nanowire. This takes into account, to lowest 
order, the hard-core repulsion of core electrons as well as the exchange energy of 
conduction electrons. When using a hard-wall confinement, the Fermi energy Ep (or 
equivalently the Fermi wavelength Xp) is the only parameter entering the NFEM. As Ep 
is material dependent and experimentally accessible, there is no adjustable parameter. 

^ Note however, that the error introduced by using hard-walls instead of a more realistic soft-wall 
confining potential can be essentially corrected for by placing the hard-wall a finite distance outside 
the wire surface, thus compensating for the over-confinement ([Garcia-Martin et al. 1996| ). 
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This pleasant feature needs to be abandoned in order to model different materials more 
realistically. Different kinds of appropriate surface boundary conditions are imaginable 
in order to model the behaviour of an incompressible fluid and to fit the surface 
properties of various metals. This will be discussed in detail in Sec. 13. 5[ 

A more refined model of a nanocontact would consider effects of scattering 
from disorder (Biirki et al. 1999a, b) and electron-electron interaction via a Hartree 
approximation f IStafford et al. 2000ai Zhang et al. 2005 ). The inclusion of disorder in 
particular leads to a better quantitative agreement with transport measurements, but 
does not change the cohesive properties qualitatively in any significant way, while 
electron-electron interactions are found to be a small correction in most cases. As a 
result, efforts to make the NFEM more realistic do not improve it significantly, while 
removing one of its main strengths, the absence of any adjustable parameters. 

The major shortcoming of the NFEM is that its applicability is limited to good 
metals having a nearly spherical Fermi surface. It is best suited for the (highly reactive) 
s-orbital alkali metals, providing a theoretical understanding of the important physics 
in nanowires. The NFEM has also been proven to qualitatively (and often semi- 
quantitatively) describe noble metal nanowires, and in particular, gold. Lately, it has 
been shown that the NFEM can even be applied (within a certain parameter range) 
to describe the multivalent metal aluminum, since Al shows an almost spherical Fermi 
surface in the extended-zone scheme. The NFEM is especially suitable to describe 
shell effects due to the conduction-band s-electrons, and the experimental observation 
of a crossover from atomic-shell to electron-shell effects with decreasing radius in both 
metal clusters ( IMartin 1996P and nanowires ( lYanson et al. 20"0T|l justifies a posteriori 
the use of the NFEM in the later regime. Naturally, the NFEM does not capture effects 
originating from the directionality of bonding, such as the effect of surface reconstruction 
observed for Au. For this reason it cannot be used to model atomic chains of Au atoms, 
which are currently extensively studied experimentally. Keeping these limitations in 
mind, the NFEM is applicable within a certain range of radius, capturing nanowires 
with only very few atoms in cross-section up to wires of several nanometers in thickness, 
depending on the material under consideration. 



3. Formalism of the NFEM 

3.1. Scattering matrix formalism 

A metal nanowire represents an open system connected to metallic electrodes at each 
end. These macroscopic electrodes act as ideal electron reservoirs in thermal equilibrium 
with a well-defined temperature and chemical potential. When treating an open system, 
the Schrodinger equation is most naturally formulated as a scattering problem. The 
basic idea of the scattering approach is to relate physical properties of the wire with 
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transmission and reflection amplitudes for electrons being injected from the leadsj^ 

The fundamental quantity describing the properties of the system is the energy- 
dependent unitary scattering matrix S{E) connecting incoming and outgoing asymptotic 
states of conduction electrons in the electrodes. For a quantum wire, S{E) can be 
decomposed into four submatrices Sap^E), a, /3 = 1,2, where 1 (2) indicates the left 
(right) lead. Each submatrix S^i^iE) determines how an incoming eigenmode of lead j3 
is scattered into a linear combination of outgoing eigenmodes of lead a. The eigenmodes 
of the leads are also referred to as scattering channels. 

The formulation of electrical transport in terms of the scattering matrix was 
developed by Landauer and Biittiker: The (linear response) electrical conductance G 
can be expressed as a function of the submatrix S21 which describes transmission from 
the source electrode 1 to the drain electrode 2 and is given by (IDatta 1995^ 

G = ^ldE ^^Tn {SI{E)S,^{E)} . (1) 

Here f{E) = {exp[l3{E — fi)] + 1}^^ is the Fermi distribution function for electrons 
in the reservoirs, /3 = (kBT)^^ is the inverse temperature and n is the electron 
chemical potential, specified by the macroscopic electrodes. The trace Tri sums over all 
eigenmodes of the source. 

The appropriate thermodynamic potential to describe the energetics of an open 
system is the grand canonical potential 

n = -l [dE D{E) ln[l + e-^^^-^'A , (2) 



(3 

where D{E) is the electronic density of states (DOS) of the nanowire. Notably, the 
DOS of an open system may also be expressed in terms of the scattering matrix as 
fIDashen et al. 19691) 



D(E).^^T.[SKE)%-%S(E)Y (3) 

where Tr sums over the states of both electrodes. This formula is also known as Wigner 
delay. Note that Eqs. ([I]), 02]), and (E]) include a factor of 2 for spin degeneracy. 

Thus, once the electronic scattering problem for the nanowire is solved, both 
transport and energetic quantities can be readily calculated. 



3.2. WKB approximation 

For an axially symmetric constriction aligned along the 2;-axis, as depicted in Fig. [H 
its geometry is characterized by the 2;-dependent radius R{z). Outside the constriction, 
the solutions of the Schrodinger equation decompose into plane waves along the wire 

^ Phase-coherence is assumed to be preserved in the wire (a good approximation given the size of the 
system compared to the inelastic mean-free-path) and inelastic scattering is restricted to the electron 
reservoirs only. 
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Figure 1. Upper- left part: Sketch of a nano-constriction. Within the adiabatic 
approximation, transverse and longitudinal motions are separable: the motion in the 
transverse direction is quantized, while in the longitudinal direction the electrons move 
in a potential created by the transverse energies [see Eq. (jB])]. Lower-left part: Sketch 
of transverse energies for different transverse channels ni, 712 , and 71-3 as a function of 
the z-coordinate. Channel rii is transmitted through the constriction as its maximum 
transverse energy is smaller than the Fermi energy, channel n2 is partly transmitted, 
and channel ns is almost totally reflected. Right part: Density plots of |5'„(r, 0)p for 
the three eigenmodes depicted on the lower-left part, corresponding to five states due 
to degeneracies of energies and En^ ■ 



and discrete eigenmodes of a circular billiard in transverse direction. The eigenenergies 
Efj^y of a circular billiard are given by 

E = — '-^ (A) 

where the quantum number 7^1, is the z/-th root of the Bessel function of order /i and 
Rq is the radius of the wire outside the constriction. In cylindrical coordinates r, 93, and 
the asymptotic scattering states read 

v[/^,(r, if, z) ~ e±*'=-^+*^'^ J^^(7^.r/i?o) , (5) 



where k^^{E) = \j2rae{E — Efj_^) /h^ is the longitudinal wavevector. In the following, 
we use multi-indices n = {^v) in order to simplify the notation. 

If the constriction is smooth, i.e. \dR/dz\ ^ 1, one may use an adiabatic 
approximation. In the adiabatic limit, the transverse motion is separable from the 
motion parallel to the 2;-axis even in the region of the constriction, and the channel index 
n of an incoming electron is preserved throughout the wire. Accordingly, Eqs. (jlj) and ([5]) 
remain valid in the region of the constriction, with R^ replaced by R{z). The channel 
energies become functions of z, En{z) = h^'^'^/2meR{zY , as is sketched in the lower 
part of Fig. [H and act as a potential barrier for the effective one-dimensional scattering 
problem in channel n. The corresponding Schrodinger equation for the longitudinal part 
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$ of the wave function reads 



ff^ 2777 

■^n{z) + ^[E-E.^{z)]^n{z) 







and is solved within the WKB approximation (see, e.g., (IMessiah 1999^ ) by 
<l>„(r, V9, z) ~ exp 



±i r K{E,z')dz' 
Jo 



(6) 



(7) 



'KiE,z) 

For a constriction of length L the transmission amplitude in channel n is then given by 
the familiar WKB barrier transmission factor 



tniE) = exp 



: / dzkn{E,z) 
Jo 



(8) 



Here Tn is the transmission coefficient of channel n and 6„ is the corresponding phase 
shift. The transmission amplitude gets exponentially damped in regions where the 
transverse energy is larger than the state total energyjj 
The full S-matrix is now found to be of the form 



S 



VTe 
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(9) 



where for simplicity of notation we have suppressed the channel indices and each of the 
entries is understood to be a diagonal matrix in the channels. Using the formulas of 

Wi we deduce 



Sec. 13.11 we may proceed to determine physical quantities. From Eq, 
that the electrical conductance at zero temperature reads, 

h 

2e^ 



G 



(10) 



5^ exp 



-2 / dze{En{z)-EF) 
Jo 



^ (Eniz) - E,) 



where the second line is obtained by using Eq. (|8]). Here 9{x) denotes the Heaviside 
step function {6{x) = 1 for a; > 0, otherwise.) The density of states is found to be 
connected with the phase shift Bn, 



D{E) 



2 dQnjE) 
vr V dE 



(11) 
(12) 

/E - E^{z) 

From the DOS, one gets the grand canonical potential in the limit of zero 
temperature as 



1 

TT 



'2me X - 



8Ei 
3A^ 



dzY.e{EF-E^{z)) 



Ejz) 
Ep 



3/2 



(13) 



which can then be used to calculate the tensile force and stability of the nanowire, as 
discussed in the following sections. 

^ This simplest WKB treatment does not correctly describe above-barrier reflection; a better 
approximation including this effect is described by Brandbyge et al. (1995) and by Glazman et al. (1988). 
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3.3. WKB approximation for non-axisymmetric wires 

The formalism presented in the previous subsection can be readily extended to non- 
axisymmetric wires. In general, the surface of the wire is given by the radius function 
r = R{ip, z), which may be decomposed into a multipole expansion 



where the sums run over positive integers. The parameterization is chosen in such a 
way that vrp(z)2 is the cross-sectional area at position z. The parameter functions Xm{z) 
and (j)rn{z) compose a vector A{z), characterizing the cross- sectional shape of the wire. 
The transverse problem at fixed longitudinal position z now takes the form 



with boundary condition Xn{R{f,z),(p;z) = for all ip G [0,27r]. This determines 
the transverse eigenenergies En{z) = En{p{z) , A{z)) which now depend on the cross- 
sectional shape through the boundary condition. With the cross-section parametrization 
f lT^ . their dependence on geometry can be written as 



where the shape- dependent functions 7ri(A) remain to be determined. In general, and 
in particular for non-integrable cross-sections, this has to be done by solving Eq. (fTSll 
numerically ( lUrban et al. 20061) . 

The adiabatic approximation (long-wavelength limit) implies the decoupling of 
transverse and longitudinal motions. One starts with the ansatz \&(r, y9,2;) = 
x(r, 2;)$(z) and neglects all ^-derivatives of the transverse wavefunction x- Again 
one is left with a series of effective one-dimensional scattering problems (Eq. [6]) for the 
longitudinal wave functions $ri(-2), in which the transverse eigenenergies En{p{z), A(z)) 
act as additional potentials for the motion along the wire. These scattering problems 
can again be solved using the WKB approximation and Eqs. ( ITTi) and ( |T3i) apply. 

3.4. Weyl- expansion 

Semiclassical approximations often give an intuitive picture of the important physics 
and, due to their simplicity, allow for a better understanding of some general features. A 
very early analysis of the density of eigenmodes of a cavity with reflecting walls goes back 
to Weyl (1911) who proposed an expression in terms of the volume and surface area of the 
cavity. His formula was later rigorously proved and further terms in the expansion were 
calculated. Quite generally, we can express any extensive thermodynamic quantity as 
the sum of such a semiclassical Weyl expansion, which depends on geometrical quantities 
such as the system volume V, surface area S, and integrated mean curvature C, as well as 




(14) 




(15) 




(16) 
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an oscillatory shell-correction due to quantum-size effects (IBrack 1997p . In particular, 
the grand-canonical potential (|2]) can be written as 

n = -ioV + (JsS --fsC + 6n, (17) 

where the energy density u, surface tension coefficient as, and curvature energy 7^ are, 
in general, material- and temperature-dependent coefficients. On the other hand, the 



shell correction Sfl can be shown, based on very general arguments (Strutinsky 1968 



Zhang et al. 2005), to be a single-particle effect, which is well described by the NFEM. 



3.5. Material dependence 

Within the NFEM there is only one parameter entering the calculation apart from the 
contact geometry: the Fermi energy Ep, which is material dependent and in general well 
known (see Tab. [1]). Nevertheless, the energy cost of a deformation due to surface and 
curvature energy, which can vary significantly for different materials, plays a crucial 
role when determining the stability of a nanowire. Obviously, when working with a 
free-electron model, contributions of correlation and exchange energy are not included, 
while they are known to play an essential role in a correct treatment of the surface 



energy (Lang 1973). Using the NFEM a priori implies the macroscopic free energy 
density uo = 2EFkp/15TT'^, the macroscopic surface energy cts = Epkp/lGn, and the 
macroscopic curvature energy 7^ = 2Epkp/9TT'^. When drawing conclusions for metals 
having surface tensions and curvature energies that are rather different from these values, 
one has to think of an appropriate way to include these material-specific properties in 
the calculation. 

A convenient way of modeling the material properties without losing the pleasant 
features of the NFEM is via the implementation of an appropriate surface boundary 
condition. Any atom-conserving deformation of the structure is subject to a constraint 
of the form 

Af = klV - 7]s k^S + r]c kpC = const. (18) 

This constraint on deformations of the nanowire interpolates between incompressibility 
and electroneutrality as side conditions, that is between volume conservation (77^ = r/c = 
0) and treating the semiclassical expectation value for the charge Qweyi (IBrack 19971) as 
an invariant {rj^ = Sir /8, r]c = 1). 

The grand canonical potential of a free-electron gas confined within a given 
geometry by hard- wall boundaries, as given by Eq. ( fT7|) . changes under a deformation 
by 

An= -cjAV + ct, A5 -7, AC + A[5fi] (19) 
= - ^AAf + {as - I^Vs)AS - (7, - ^Vc)AC + A[5n], 

where the constraint (1181) was used to eliminate V. Now the prefactors of the change 
in surface AS and the change in integrated mean curvature AC can be identified as 
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Table 1. Material parameters (jAshcroft and Mermin 19761 IPerdew et al. 199T|l of 
several monovalent metals: Fermi energy Ep^ Fermi wavevector kp, surface tension 
CTs, and curvature energy 7^, along with the corresponding values of 77^ and rjc- The 
last column gives the corresponding values for the multivalent metal Al (see discussion 
in Sec. Ell) Adapted from (jUrban et al. 2006^ . 

effective surface tension and curvature energy, respectively. Tliey can be adjusted to 
fit a specific material's properties by an appropriate choice of tfie parameters r]s and rjc 
(see Tab. [I]). 

4. Conductance and Force 

The formalism presented in the previous section can now be applied to a specific wire 
geometry (IStafford et al. 19W|) . namely a cosine constriction, 

R{z) = + cos {-J-) , (20) 

of a cylindrical wire. One is interested in the mechanical properties of this metallic 
nanoconstriction in the regime of conductance quantization. The necessary condition to 
have well-defined conductance plateaus in a three-dimensional constriction was shown 
fITorres et al. 1994p to be [dR/dzf < 1. In this limit, one may employ the adiabatic 
and WKB approximations and evaluate the expressions obtained in Sec. 13.21 

4.1. Conductance 

The conductance is obtained from Eq. (fTOj) . As the transmission amplitudes Tn vary 
exponentially from 1 to when the transverse energy of the respective channel at the 
neck of the constriction traverses the Fermi energy, this results in a steplike behavior 
of the conductance with almost fiat plateaus in between. This is the phenomenon 
of conductance quantization, which is observable even at room temperature for noble 
metal nanowires due to the large spacing of transverse energies (of order 1 eV for Au, 
to compare to ksT 10^'^eV at room temperature). The upper panel of Fig. Et^a) 
shows the conductance obtained with an improved variant of the WKB approximation 
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Figure 2. Electrical conductance G and tensile force F as function of elongation of a 
nanowire, (a) calculated in a WKB approximation for a cosine constriction (adapted 
from Stafford et al. 1997), and (b) measured by Rubio et al. (1996) in an experiment 
on gold (data courtesy of N. Agrait). 



( IBrandbyge et al. IQQSj [Gl azman et al. 19881) for the geometry ( l20l) . The conductance 



as a function of elongation shows the expected stephke structure and the step heights 
are 2e^/h and integer multiples thereof (the multiplicity depends on the degeneracy of 
the transverse modes). An ideal plastic deformation was assumed, i.e. the volume of 
the constriction was held constant during elongation]^ 

4.2. Force 

If the wire elongation is slow enough, the electron gas has time to adjust to the wire 
shape changes during the deformation, and is thus always in equilibrium]^ Under these 
conditions, the tensile force can be computed from the grand canonical potential, given 
in the WKB approximation by Eq. (fT3|) . as F = —dfl/dL\^. 

The lower panel of Fig. |2|^a) shows the tensile force for the cosine constriction ( l20l) . 
The correlations between the force and conductance are striking: \F\ increases along 
the conductance plateaus, and decreases sharply when the conductance drops. The 
constriction becomes unstable when the last conductance channel is cut off, which is 

Note the different abscissae for the theoretical and experimental graphs in Fig. [21 While both 
contact diameter and elongation are a measure of the deformation of the contact, the latter is more 
easily accessible to experiments. The contact diameter on the other hand is the natural independent 
geometric variable which is used for the theoretical graph, since its relation to elongation is very model- 
dependent. 

^ The experimental elongation speed is of order Inm/s, which is compatible with this assumption. 
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indicated by a positive tensile force. Some transverse channels are quite closely spaced, 
and in these cases, the individual conductance plateaus [e.g., G/Gq = 14, 15, 19, 21] 
and force oscillations are difficult to resolve. The force oscillations are found to have 
an amplitude Ep/\p (i.e. ~ 1.7 nN for gold, consistent with experimental observations) 
independent of the chosen geometry (circular and quadratic wires and cosine and 
parabolic constrictions were tested), and to persist to very large conductances. 

These results can be understood within the intuitive picture of a conductance 
channel as a delocalized metallic bond. The increase of \F\ along the conductance 
plateaus and the rapid decrease at the conductance steps can then be interpreted as 
stretching and breaking of these bonds. Note that within the NFEM, the correlation 
between conductance changes and force oscillations comes from a pure quantum-size 
effect and not from atomic rearrangements. 

The comparison of theoretical predictions with experimental data by Rubio et 
al. (1996), plotted in Fig. |2]^b), shows very good qualitative agreement and underlines 
that, although the NFEM is simple, it already gives a good qualitative description of 
the experimental findings. Extensions of the NFEM including structural dynamics of 
the wire (Biirki et al. 2003, 2007a,b), which are able to calculate the shape of the wire 
at all steps of elongation, find that instabilities accelerate the conductance and force 
jumps, making the theory even more similar to the experiment. 

In addition, it is easy to show that in the NFEM, the tensile force is invariant under 
a stretching of the geometry i?(2;) R{Xz) , so that F = {ef/^f) /(AL/Lq, /ci? -R), where 
f{x,y) is a dimensionless function, i.e. the force oscillations are universal, and thus do 
not depend on the details of the wire shape, or precisely how it deforms. Nonuniversal 
corrections to F occur in very short constrictions, for which the adiabatic approximation 
breaks down. 

5. Linear Stability Analysis 

Metal nanowires are of great interest for nanotechnology, since they may serve as 
conductors in future nanocircuits. In particular, one would like to know whether a 
nanowire of given length and radius remains stable at a given temperature. 

At first sight, an atomistic approach seems to be more "realistic" than the NFEM 
and well suited to answer this question. But molecular dynamics (MD) simulations 
conceptually are not able to avoid the surface-tension-driven Rayleigh instability of long 
nanowires. Since quantum-size effects from the electron confinement are not properly 
taken into account, MD simulations fail to give an explanation for the electronic shell 
and supershell effects. On the other hand, atomistic quantum calculations using, e.g., 
the local-density approximation, are restricted to such small systems that their results 
can not really be disentangled from finite-size effects (IStafford et al. 2000bp . Therefore, 
to date, a stability analysis within the framework of the NFEM (and generalizations 
thereof) is the only approach able to correctly include the effects of electron-shell filling 
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and thereby shed hght on the puzzhng stabihty of long metal nanowires. 

The geometry of a wire of uniform cross section aligned along the z-axis is 
characterized by the cross-sectional area A = vrp^, and a set of dimensionless parameters 
determining the shape, which compose a vector A (cf. Eq. [T^ . A small 2;-dependent 
perturbation of a wire of length L and initial cross section (p, A) can be written in terms 
of a Fourier series as 

p{z) =p + e6p{z) = p + 6j2pge"", 

A{z) = A + 6 6A{z) = A + eY.^qe'"", (21) 

where the dimensionless small parameter e sets the size of the perturbation|§ 

The energetic cost of a small deformation of the wire can be calculated by expanding 
the grand canonical potential as a series in the parameter e, 

n = + e + n^'^ + oie^). (22) 

A nanowire with initial cross-section (p, A) is energetically stable at temperature T 
if and only if fiW(p, A,T) = and n^'^\p,A,T) > for every possible deformation 
{6p,6A) satisfying the constraint ( ITSjl . 



5.1. Rayleigh- Instability 

It is instructive to forget about quantum-size effects for a moment and to perform a 
stability analysis in the classical limit. For simplicity, one can restrict oneself to axial 
symmetry (i.e. A = 0). In the classical limit the grand canonical potential is given by 
the leading order terms of the Weyl approximation, flwoyi = V + asS, and changes 
under the perturbation (12T!) by 

- 2tt {pu - as) Po£ + T^^\Pq\'^ -u + q^pasE^. (23) 



'Wcyl _ 



Because of the constraint (fT8|) on possible deformations, po can be expressed in terms of 
the other Fourier coefficients. Volume conservation, e.g., implies po = — (e/2p) Z^g^^o |PgP 
and 

^ ^ ^ |p,r(p'g' - l)e^ (24) 

which has to be positive in order to ensure stability. Since q is restricted to integer 
multiples of 27r/L, stability requires L < 27rp. This is just the criterion of the classical 
Rayleigh instability ( IChandrasekhar 19811) : A wire longer than its circumference is 
unstable and likely to break up into clusters due to surface tension. 

^ Assuming periodic boundary conditions, the perturbation wave vectors q must be integer multiples 
of 271 /L. In order to ensure that p{z) and A(z) are real, we have p^q = and = A*. 
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5.2. Quantum-mechanical stability analysis 



The crucial ingredient to the stabihzation of metal nanowires is the oscillatory shell 
correction 6Q to the grand canonical potential ( 1T71) which is due to quantum-size effects. 
This shell correction can be accounted for by a quantum-mechanical stability analysis 
based on the WKB-approximation introduced in Sec. 13.21 The use of this approximation 
can be justified a posteriori by a full quantum calculation (Urban et al. 2003, 2007) 
which shows that the structural stability of metal nanowires is indeed governed by 
their response to long-wavelength perturbations. The response to short-wavelength 
perturbations on the other hand controls a Peierls-type instability characterized by the 
opening of a gap in the electronic energy dispersion relation. This quantum mechanical 
instability, which is missing in the semiclassical WKB approximation, in fact limits the 
maximal length of stable nanowires. Nevertheless, if the wires are short enough, and/or 
the temperature is not too low, the full quantum calculation essentially confirms the 
semiclassical results. 

A systematic expansion of Eq. ( fT3ll yields 



(25) 




L/Xf ^ 
L/Xf ~~ 

where the elements of the matrix A in Eq. fl26|) are given by 
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(27) 



Ep — En 



Here E'^ denotes the gradient of En with respect to A and E'^^ is the matrix of second 
derivatives. The bar indicates evaluation at (p. A). 

The number of independent Fourier coefficients in Eq. f l2T|) is restricted through 
the constraint ( fTSl) on allowed deformations. Hence, after evaluating the change 
of the geometric quantities V, 5, and C due to the deformation, we can use Eq. 
( ITSjl to express po in terms of the other Fourier coefficients, yielding an expansion 
Po = Po°'* +^ Po^'* + ^(^^)- This expansion then needs to be inserted in Eqs. ( l25l) and ( l26l) . 
thereby modifying the first-order change of the energy QS^'> and the stability matrix A 
(lUrban et al. 2006|1 . 
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Figure 3. WKB stability coefficient, calculated using a constant-volume constraint. 
The sharp negative peaks at the opening of new channels (i.e. when kpp = 7„) are 
smeared out with increasing temperature T. 



Stability requires that the resulting modified stability matrix A be positive definite. 
Results at finite temperature are obtained essentially in a similar fashion, by integrating 
Eq. ([2]) numerically. 

5.3. Axial symmetry 

A straightforward application of the method outlined above is the stability analysis of 
cylindrical wires with respect to axisymmetric volume conserving perturbations. In this 
specific case A(z) = and po = ~(^/2p) J2q^o IPgP- Therefore, Eq. ( I25l) takes the form 
fi(^) = and Eq. ( 126|) simplifies to read 

= EfY.\pJp\ (^ip), (28) 



where the stability coefficient a{p) = App reads (lUrban et al. 2006|) 

47^ 



4. 



>y2 

2 'n 



(29) 



Axial symmetry implies the use of the transverse eigenenergies En/Ep = {'fn/kpf))'^, 
cf. Eq. @. This result, valid for zero temperature, is plotted as a function of radius 
in Fig. |3] together with a numerical result at finite temperature. Sharp negative peaks 
at the subband thresholds, i.e. when pkp = 7„, indicate strong instabilities whenever 
a new channel opens. On the other hand, a is positive in the regions between these 
thresholds giving rise to intervals of stability that decrease with increasing temperature. 
These islands of stability can be identified with the "magic radii" found in experiments. 
As will be shown below, one has to go beyond axial symmetry in order to give a full 
explanation of the observed conductance histograms of metal nanowires. 
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G / Gq a, 



A 2 T^.^^/Tp 



2 1.72 

5 1.33 

9 1.22 

29 1.13 

59 1.11 

72 1.08 

117 1.06 

172 1.06 



0.26 0.50 

0.14 0.49 

0.10 0.50 

0.06 0.54 

0.05 0.49 

0.04 0.39 

0.03 0.55 

0.03 0.50 



Table 2. Most stable deformed wires with quadrupolar cross sections. The first 
column gives the quantized conductance of the corresponding wire. Both the aspect 
ratio a and the value of the deformation parameter A2 are given. The maximum 
temperature of stability T^-^ax is given for each wire. In all cases the surface tension 
was set to 0.22 N/m, corresponding to Na. Adapted from (jUrban et al. 20061) . 



5.4- Breaking axial symmetry 

It is well known in the physics of crystals and molecules that a Jahn- Teller deformation 
breaking the symmetry of the system can be energetically favorable. In metal clusters, 
Jahn- Teller deformations are also very common, and most of the observed structures 
show a broken spherical symmetry. By analogy, it is natural to assume that for 
nanowires, too, a breaking of axial symmetry can be energetically favorable, and lead 
to more stable deformed geometries. 

Canonical candidates for such stable non-axisymmetric wires are wires with a 
cos(m(y9)-deformed cross section (i.e. having m-fold symmetry), a special case of Eq. f|T^ 
with only one non-zero A^. The quadrupolar deformation (m = 2) is expected to be 
the energetically most favorable of the multipole deformation^ since deformations with 
m > 2 become increasingly costly with increasing m, their surface energy scaling as m^. 

The results of Sec. 15.21 can straightforwardly be used to determine stable 
quadrupolar configurations by intersection of the stationary curves, il^^\p,\2)\j\f = 0, 
and the convex regions, 

n^^\p,X2)W > 0. The result is a so-called stability diagram 
which shows the stable geometries (at a given temperature) in configuration space, that 
is a as function of the geometric parameters p and A2. An example of such a stability 
diagram is shown later in Fig. [6] for the case of aluminum, discussed in Sec. 15.71 Results 
for all temperatures can then be combined, thus adding a third axis (i.e. temperature) 
to the stability diagram. Finally, the most stable configurations can be extracted, 
defined as those geometries that persist up to the highest temperature compared to 
their neighboring configurations. 

The dipole deformation (m — 1) corresponds, in leading order, to a simple translation, plus higher- 
order multipole deformations. Therefore the analysis can be restricted to m > 1. 
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Table [2]lists the most stable deformed sodium wires with quadrupolar cross section, 
obtained by the procedure described above. The deformation of the stable structures is 
characterized by the parameter A2 or equivalently by the aspect ratio 

a = ^, ' -. (30) 

V^l - Ai/2 - A2 

Clearly, nanowires with highly-deformed cross sections are only stable at small 
conductance. The maximum temperature up to which the wires remain stable, given 
in the last column of Tab. 121 is expressed in units of Tp := Tp/{kFp)- The use of this 
characteristic temperature reflects the temperature dependence of the shell correction 
to the wire energy ( lUrban et al. 2006]) . 

Deformations with higher m cost more and more surface energy. Compared to 
the quadrupolar wires, the number of stable configurations with 3-, 4-, 5-, and 6-fold 
symmetry, their maximum temperature of stability, and their size of the deformations 
involved, all decrease rapidly with increasing order m of the deformation. For m > 6 
no stable geometries are known. All this reflects the increase in surface energy with 
increasing order m of the deformation. 



5. 5. General stability of cylinders 

It is possible to derive the complete stability diagram for cylinders, i.e., to determine the 
radii of cylindrical wires that are stable with respect to arbitrary small, long-wavelength 
deformations (lUrban et al. 20061) . ^^st sight, considering arbitrary deformations. 



and therefore theoretically an infinite number of perturbation parameters seems a 
formidable task. Fortunately, the stability matrix A for cylinders is found to be 
diagonal, and therefore the different Fourier contributions of the deformation decouple. 
This simplifies the problem considerably, since it allows to determine the stability of 
cylindrical wires with respect to arbitrary deformations through the study of a set of 
pure m-deformations, i.e. deformations as given by Eq. ( IT4l) with only one non-zero A^. 

Figure m shows the stable cylindrical wires in dark gray as a function of temperature. 
The surface tension was fixed at the value for Na, see Tab. [H The stability diagram 
was obtained by intersecting a set of individual stability diagrams allowing cos(my9)- 
deformations with m < 6. This analysis confirms the extraordinary stability of a set of 
wires with so called "magic radii". They exhibit conductance values G/Gq = 1, 3, 6, 12, 
17, 23, 34, 42, 51,... It is noteworthy that some wires that are stable at low temperatures 
when considering only axisymmetric perturbations, e.g., G/Gq = 5, 10, 14, are found to 
be unstable when allowing more general, symmetry-breaking deformations. 

The heights of the dominant stability peaks in Fig. HI exhibit a periodic modulation, 
with minima occurring near G/Gq = 9, 29, 59, 117, ... The positions of these minima are 
in perfect agreement with the observed supershell structure in conductance histograms 
of alkali metal nanowires OYanson et al. 20"00ll . Interestingly, the nodes of the supershell 
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Figure 4. Stability of metal nanocylinders versus electrical conductance and 
temperature. Dark gray areas indicate stability with respect to arbitrary small 
deformations. Temperature is displayed in units of Tp = Tp/kpp (see text). 
The surface tension was taken as 0.22 N/m, corresponding to Na. Adapted from 
(jUrban et al. 2006^ . 



structure, where the shell effect for a cylinder is suppressed, are precisely where the 
most stable deformed nanowires are predicted to occur (see discussion above). Thus 
symmetry breaking distortions and the supershell effect are inextricably linked. 

Linear stability is a necessary — but not a sufficient — condition for a nanostructure 
to be observed experimentally. The linearly stable nanocylinders revealed in the above 
analysis are in fact metastahle structures, and an analysis of their lifetime has been 
carried out within an axisymmetric stochastic field theory by Biirki et al. (2005a). There 
is a strong correlation between the height of the stable fingers in the linear stability 
analysis and the size of the activation barriers Aii^, which determines the nanowire 
lifetime r through the Kramers formula r = tq exp(Aii^/fcBT). This suggests that the 
linear stability analysis, with temperature expressed in units of Tp = Tp/kFP, provides 
a good measure of the total stability of metal nanowires. In particular, the "universal" 
stability of the most stable cylinders is reproduced, wherein the absolute stability of 
the magic cylinders is essentially independent of radius (aside from the small supershell 
oscillations) . 

5.6. Comparison with experiments 

A detailed comparison between the theoretically most stable structures and experimental 
data for sodium is provided in Fig. |5l For each stable finger in the linear stability 
analysis its mean conductance is extracted and plotted as a function of its index number, 
together with experimental data by Yanson et al. (1999). This comparison shows that 
there is a one-to-one relation between observed conductance peaks and theoretically 
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Figure 5. Comparison of tlie experimental shell structure for Na, taken from Yanson 
et al. (1999), with the theoretical predictions of the most stable Na nanowires. Non- 
axisymmetric wires are labeled with the corresponding aspect ratio a. Adapted from 
(|Urban et al. 2004ap . 



stable geometries which in particular allows for a prediction of the cross-sectional shape 
of the wires. This striking fit is only possible when including non-axisymmetric wires, 
which represent roughly 25% of the most stable structures and which are labeled by 
the corresponding aspect ratios a in Fig. [5l The remaining 75% of the principal 
structures correspond to the magic cylinders. The role of symmetry in the stability 
of metal nanowires is thus fundamentally different from the case of atomic nuclei or 
metal clusters, where the vast majority of stable structures have broken symmetry. The 
crucial difference between the stability of metal nanowires and metal clusters is not the 
shell effect, which is similar in both cases, but rather the surface energy, which favors 
the sphere, but abhors the cylinder. 

Besides the geometries entering the comparison above, the stability analysis also 
reveals two highly deformed quadrupolar nanowires with conductance values of 2 Gq 
and 5 Gq, cf. Tab. El They are expected to appear more rarely due to their reduced 
stability relative to the neighboring peaks, and their large aspect ratio a that renders 
them rather isolated in configuration space|§ Nevertheless they can be identified by a 
detailed analysis of conductance histograms of the alkali metals (lUrban et al. 2004bll . 

^ A nanowire produced by pulling apart an axisymmetric contact has a smaller probability to transform 
into a highly deformed configuration than into a neighboring cylindrical configuration. 
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Figure 6. Stability diagram for Al wires at fixed temperature T — 0.45 Tp. Thick 
lines mark stable wires in the configuration space of rms radius p and deformation 
parameter A2. The dashed box emphasizes a series of very stable superdeformed wires, 
whose peanut-shaped cross section is shown as an inset. This sequence was recently 
identified experimentally (jMares tt al. 2007"|) . 



5.7. Material dependence 

Results for different metals are similar in respect to the number of stable configurations 
and the conductance of the wires. On the other hand, the deviations from axial 
symmetry and the relative stability of Jahn- Teller deformed wires is sensitive to the 
material-specific surface tension and Fermi temperature. The relative stability of the 
highly deformed wires decreases with increasing surface tension as/ {Epk^), measured 
in intrinsic units, and this decrease becomes stronger with increasing order m of the 
deformation. Therefore, for the simple s-orbital metals under consideration (Tab. [1]), 
deformed Li wires have the highest, and Au wires have the lowest relative stability 
compared to cylinders of "magic radii. "§ Notable in this respect is Aluminum with 
as = O.OOnEpk'j^, some five times smaller than the value for Au. Aluminum is a trivalent 
metal, but the Fermi surface of bulk Al resembles a free-electron Fermi sphere in the 
extended-zone scheme. This suggests the applicability of the NFEM to Al nanowires, 
although the continuum approximation is more severe than for monovalent metals. 

Recent experiments ( IMares et al. 2007]) have found evidence for the fact that the 



stability of aluminum nanowires also is governed by shell filling effects. Two magic 
series of stable structures have been observed with a crossover at G ^ 406*0 and the 
exceptionally 

^ Concerning the absolute stability, we have to consider that the lifetime of a metastable nanowire also 
depends on the surface tension IjBiirki et al. 2005a|) . 
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stable structures have been related to electronic and atomic shell effects, 
respectively. Concerning the former, the NFEM can quantitatively explain the 
conductance and geometry of the stable structures for wires with G > 12Go and 
there is a perfect one-to-one correspondence of the predicted stable Al nanowires and 
the experimental electron-shell structure. Moreover, an experimentally observed third 
sequence of stable structures with conductance G/Gq ~ 5,14,22 provides intriguing 
evidence for the existence of "super deformed" nanowires whose cross sections have an 
aspect ratio near 2:1. Theoretically, these wires are quite stable compared to other highly 
deformed structures and, more importantly, are very isolated in configuration space, 
as illustrated in the stability diagram shown in Fig. |6l This favors their experimental 
detection if the initial structure of the nanocontact formed in the break junction is rather 
planar with a large aspect ratio since then it is likely that the aspect ratio is maintained 
as the wire necks down elastically. Aluminum is unique in this respect and evidence of 
superdeformation has not been reported in any of the previous experiments on alkali 
and noble metals, presumably because highly-deformed structures are intrinsically less 
stable than nearly axisymmetric structures, due to their larger surface energy. 

6. Summary and discussion 

In this chapter we have given an overview on the Nanoscale Free-Electron Model, 
treating a metal nanowire as a non-interacting electron gas confined to a given geometry 
by hard- wall boundary conditions. At first sight, the NFEM seems to be an overly simple 
model, but closer study reveals that it contains very rich and complex features. Since its 
first introduction in 1997, it has repeatedly shown that it captures the important physics 
and is able to explain qualitatively, when not quantitatively, many of the experimentally 
observed properties of alkali and noble metal nanowires. Its strengths compared to other 
approaches are, in particular, the absence of any free parameters and the treatment of 
electrical and mechanical properties on an equal footing. Moreover, the advantage of 
obtaining analytical results allows the possibility to gain some detailed understanding 
of the underlying mechanisms governing the stability and structural dynamics of metal 
nanowires. 

The NFEM correctly describes electronic quantum size effects, which play an 
essential role in the stability of nanowires. A linear stability analysis shows that the 
classical Rayleigh instability of a long wire under surface tension can be completely 
suppressed by electronic shell effects, leading to a sequence of certain stable "magic" wire 
geometries. The derived sequence of stable cylindrical and quadrupolar wires explains 
the experimentally observed shell and supershell structures for the alkali and noble 
metals as well as for aluminum. The most stable wires with broken axial symmetry are 
found at the nodes of the supershell structure, indicating that the Jahn- Teller distortions 
and the supershell effect are inextricably linked. In addition, a series of superdeformed 
aluminum nanowires with an aspect ratio near 2:1 is found which has lately been 
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identified experimentally. A more elaborate fully quantum mechanical analysis within 
the NFEM reveals an interplay between the Rayleigh and a Peierls-type instability. The 
latter is length-dependent and limits the maximal length of stable nanowires but other 
than that confirms the results obtained by the long wavelength expansion discussed 
above. Remarkably, certain gold nanowires are predicted to remain stable even at 
room temperature up to a maximal length in the micrometer range, sufficient for future 
nanotechnological applications. 

The NFEM can be expanded by including the structural dynamics of the wire in 
terms of a continuum model of the surface diffusion of the ions. Furthermore, defects 
and structural fluctuations may also be accounted for. These extensions improve the 
agreement with experiments but do not alter the main conclusions. However, the NFEM 
does not address the discrete atomic structure of metal nanowires. With increasing 
thickness of the wire the effects of surface tension decrease and there is a crossover from 
plastic flow of ions to crystalline order, the latter implying atomic shell effects observed 
for thicker nanowires. Therefore, the NFEM applies to a window of conductance values 
between a few Go ^-nd about IOOGq) depending on the material under consideration. 

Promising extensions of the NFEM in view of current research activities are 
directed, e.g., towards the study of metal nanowires in nanoelectromechanical systems 
(NEMS) which couple nanoscale mechanical resonators to electronic devices of similar 
dimensions. The NFEM is ideally suited for the investigation of such systems since 
it naturally comprises electrical as well as mechanical properties. It is hoped that 
the generic behaviour of metal nanostructures elucidated by the NFEM can guide 
the exploration of more elaborate, material-speciflc models, in the same way that the 
free-electron model provides an important theoretical reference point from which to 
understand the complex properties of real bulk metals. 
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